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Abstract 

This paper presents an inverse problem for the nonlinear 1-d Kuramoto-Sivashinsky 
(K-S) equation. More precisely, we study the nonlinear inverse problem of retrieving the 
anti-diffusion coefficient from the measurements of the solution on a part of the boundary 
and at some positive time everywhere. Uniqueness and Lipschitz stability for this inverse 
problem are proven with the Bukhgeim-Klibanov method. The proof is based on a global 
Carleman estimate for the linearized K-S equation. 
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1 Introduction 

This paper focuses on an inverse problem that consists in the determination of a coefficient 
of a partial differential equation from the partial knowledge of a given single solution of the 
equation. For the solution of this class of problems (single-measurement inverse problems) , 
the Bukhgeim-Klibanov method was introduced in [6] (see also |121|I3] ). This method, which 
is based on Carleman estimates, allows to prove uniqueness, i.e. that each measurement 
corresponds to only one coefficient. Regarding the continuity of this inverse problem, the 
first Lipschitz stability result for a multidimensional wave equation was obtained by Puel 
and Yamamoto |14] by using a method based on |B]. Since then, this method has been 
applied to other inverse problems including Lipschitz stability for the Schrodinger equation 
[I] and Logarithmic stability for the wave equation [21 [3] . 
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This approach was extended to a parabolic equation in QJ5]. Since then, this type of 
inverse problems for parabolic equations has received a large amount of attention. The 
primary difference with respect to hyperbolic inverse problems is that parabolic problems 
are not time-reversible: therefore, an additional measurement must be added if that method 
is applied. As one can read in the discussion of the introduction of [ID], the knowledge of the 
full-state of the solution for some positive time is required . To prove the Lipschitz stability 
without this assumption, which is usually needed when global Carleman inequalities are 
used, is still an open problem. 

Recent results regarding linear parabolic problems can be found in [3] (discontinuous 
coefficient), [8] (systems), [11] (network) and the references therein. In 5, 9, 15., nonlinear 
parabolic equations were even considered. 

In this paper, the Kuramoto-Sivashinsky (K-S) equation is considered, which is a ID 
nonlinear fourth-order parabolic equation. This equation is used to model the physical 
phenomena of plane flame propagation: it describes the combined influence of diffusion and 
thermal conduction of gas on the stability of a plane flame front. In this nonlinear partial 
differential equation, the fourth-order term models the diffusion, and the second-order term 
models the incipient instabilities. To the knowledge of the authors there are no results 
in the literature concerning the determination of coefficients for this nonlinear equation. 
However, a Carleman estimate has been used to obtain the null-controllability of the K-S 
equation in reference [7] for the constant coefficient case. We consider the inverse problem 
of retrieving the anti-diffusion coefficient 7 from boundary measurements of the solution. 
Since the linearized equation is parabolic, boundary measurements are not sufficient and 
we must consider an additional measurement of the full solution for a given time To (as in 
[U [10] among others) . 

The K-S equation with non-constant coefficients describing the diffusion a = cr(x), and 
the anti-diffusion 7 = y(x), is given as 



where Q := (0, T) x (0, 1), a : [0, 1] — > R+, and the functions j/o, g, hj are the initial condition, 
the source term and the boundary data respectively. All of these terms are assumed to be 
known and compatible. 

The first result involves the local well-posedness of the nonlinear equation {TJ. A less 
regular framework can be used for this equation. However, the method applied in this paper 
requires the solution and its time-derivative to be at least in L 2 (0, T; H 4 (0, 1)). Therefore, 



< 



yt + (a(x)y xx ) xx + j(x)y xx + yy x = g, 
y(t,0) =hi(t), y(t,l) = h 2 (t), 
y»(*,0)=As(t), y x {t,l) = h 4 (t), 
2/(0, as) = yo(x), 



V(t,a:) 6 Q, 
Vie (0,T), 
Vie (0,T), 
VxG (0,1), 



(1) 



2 



let us introduce the following notations for the spaces appearing in this paper: 



y k := C([0, T];H k (0, 1)) n L 2 (0, T; H k+2 (0, 1)), for k G N; 
T := {/ G L 2 (0, T; ff 4 (0, 1))/ f t G L 2 (0, T; L 2 (0, 1))}; 

z~{zey e / Zt ey 2 }. 



(2) 



Theorem 1.1 Let 7 G H 4 (0, 1) and a G ff 4 (0, 1) be such that 



Vx G (0,1), a(x) > ao > 0. 



(3) 



Let 2/0 S H 6 (0, 1), j 6 J, and /ij G H 2 (0,T) for j — 1,...,4. Assume also that yo and 
(ft.j)j=i,...,4 satisfy the compatibility conditions j/o(0) = /ii(0), yo,x(0) = /13(C)), 2/o(l) = /i2(0), 
2/o,x(l) = ft 4 (0). 

Therefore, there exists e > stic/i t/iat i/ 



i/ien </ie ./('-S' equation JTJ /ios a unique solution y £ Z . 

Once the existence of solutions to the K-S equation has been established (see Section f2|, 
the following inverse problem is addressed: 

Is it possible to retrieve the anti-diffusion coefficient 7 = *y(x) from the measure- 
ment of y X x(t, 0) and y xxx {t,Q) on (0,T) and from the measurement of y(To,x) 
on (0, 1), where y is the solution to Equation |T} and To G (0, T)? 

A local answer for this nonlinear inverse problem is given (see section f?]). To be more spe- 
cific, let 7 be fixed. We denote by y the solution to Equation |T} with 7 replaced by 7. This 
paper focuses on the following two questions. 

Uniqueness: Do the equalities of the measurements y X x(t,0) ~ y X x(t,0) and y X xx(t, 0) = 
yxxx(t,0) for t G (0,T) and y(T ,:r) = y{T ,x) for a: G (0, 1) imply 7 = 7011 (0, 1)? 

Stability: Is it possible to estimate H7 — 7||l 2 (o,i) by suitable norms ||y(Tb,:r) — y(To,x)|| 
in space and \\y~xx(t,0) - yxx(t,0)\\, \\y X xx(t, 0) - y X xx{t, 0)|| in time? 

To answer these questions, we use the Bukhgeim-Klibanov method. First, a global Car- 
leman estimate for the linearized K-S equation with non-constant coefficients is obtained. It 
is then used to prove the primary result which can be stated as follows. 

Theorem 1.2 Let us consider a G H 4 (0, 1), 7 G H 4 (0, 1) and compatible data yo, g and hj 
regular enough such that the solutions of <[TJ) belong to i? 1 (0, T; H 4 (0, 1)). 

Let us denote y the solution of equation J7]) associated to 7 and y the solution associated 
to a fixed coefficient 7 G H 4 (0,L). We assume that there exists r > and To G (0, T) such 
that 



1 1 2/o 1 1 h<5 (0,1) < e> ||fllk<e> HMh 2 (o,t) < e for j = h ••• »4, 



(4) 



mi{\y X x(To,x)\,x€ (0, 1)} > r, 



(5) 
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Let U — {7 G H 4 (0, 1) / ||7|| loo (0,1) < Mi, ||i/||iri(o,T:H 4 (o,i)) < A^a}- T/ien f/iere exists a 
positive constant C depending on the parameters (T, M\,M 2 ,r), such that for every 7 G W, 

^ll7-7lli 2 (o,i) < \\vxx{-,Q) - yxx{-,o)\\ 2 H i {0]T) + \\vxxx(-,Q) - yxxx(-,o)\\ 2 H i(o,T) 
+ \\y( T o, •) - y( T o>-)llff4(o,i) + \\vi T o,-) - y(T ,-)\\ 

< C(\\y - 2/||Hl(0,TiH<i(0,l)) + ||y - y||ioo(0,T;Hl(0,l)))- ( 6 ) 

This two-sided inequality gathers two complementary informations, namely, the stability 
of the inverse problem (with the first estimate) and the regularity of the measurements. 
Indeed, the second estimate in ([6]) indicates that the required measurements are finite if y 
and y belong to the space H 1 (0, T; H 4 (0, 1)) and this is true if y and y are solutions provided 
by Theorem 1 1 . 1 1 under the hypothesis ((3| and (Q. 

Remark 1.3 As stated above, an internal measurement at t = To is required if this method 
is used to solve this type of problem for parabolic equations. Nevertheless, this is probably a 
technical point and there is no counter-example that demonstrates whether this assumption 
is required. In .15,, uniqueness (but not stability) is proven using a very different technique 
in an inverse problem for a parabolic equation and without any internal measurements in the 
entire space domain. 

Remark 1.4 In this paper, the boundary measurements are located at x = 0, but the result 
would be the same if we measure at x — 1 instead. Indeed, the choice of a suitable weight 
function in the proof of the Carleman estimate in Section^ is critical to impose the side of 
measurement. 

This article is organized as follows. The well-posedness result stated in Theorem 11.11 is 
proved in Section [21 A global Carleman estimate for a general K-S equation is given and 
proved in Section [3] Finally, Section [4] conatins the use of the Bukhgeim-Klibanov method 
to prove the Lipschitz stability of the inverse problem stated in Theoren ll.2l 



2 On the Cauchy problem for KS equation 

This section presents a proof of Theorem ll.ll in a more general case including time dependent 
lower-order coefficients. We consider the following K-S system 



y t + {ar{x)y X x)xx + r ){x)y xx + G\y x + G 2 y + yy x = g, V(t, x) e Q, 

y(t,0) = hi(t), y(t,l) =Aa(t), Vte(0,T). 

y x (t,0) = h 3 (t), y x (t,l) = h i (t), Vte(0,T), 

y(0,x)=y (x), Vx 6 (0, 1), 



(7) 



where Gi,G 2 belongs to J/ 1 (0, T; H 4 (0, 1)), g G T and y Q G ff 6 (0, 1) is compatible with 
hj G H 2 (0, T) for j = 1, ... , 4. Recall that the coefficients satisfy 7 G # 4 (0, 1), a G H 4 (0, 1) 



4 



and hypothesis (J3J). 



First, the main part of the linear differential operator is utilized in the next proposition. 
Proposition 2.1 Let zo G H G n Hg(0, 1) and f 6 T . Then, the following equation 
zt + (ct(x)z X x)xx = f, V(t,z)eQ, 



(8) 



«(t,0)=0, «(t,l)=0, Vt£(0,T), 
z x (t,0) = 0, z x (t, 1) = 0, Vte(0,T), 
2(0, a;) =«o(aj), Vie (0,1), 

/ias a unique solution z G Z and there exists C > suc/i £/iat 

INU<c (||/||^ + w„.). 

Proof. The operator 

H 4 n H 2 (0, 1) C L 2 (0, 1) — > L 2 (0, 1) 

tu i— > (<r(a;)w"(a;))", 
is simultaneously positive, coercive and self-adjoint. Moreover, its inverse is compact: thus, it 
generates a strongly continuous semigroup in L 2 (0, 1). Therefore, for each zo G H 4 (~)Hq(0, 1) 
and / G C 1 ([0, T]; L 2 (0, 1)), Equation © has a unique solution 2 G C([0, T]; H 4 flHi (0, l))n 
C l (t0,ri;L a (0,l)). 

We will demonstrate that the solutions z £ Z (refer to the notation introduced in ©), 
can be obtained by taking zq and / sufficiently regular. 

We now search for some energy estimates that indicate the space that the solutions lie on 
depending on the regularity of the data. Suppose that there are solutions sufficiently regular 
to perform the following computations. Equation © is multiplied by z and integrated over 
(0, 1). Some integrations by parts give 



d 
~dt 



z\ 2 )+ I \z xx \ 2 <Cl f\.f\ 2 + C\z\ 2 ). (!)) 
o / Jo \Jo Jo 



Throughout this paper, C denotes a positive constant that varies from line to line. Using 
GronwalPs lemma, we obtain 



i / r r r 1 

2 - s < I II I j? 1 2 | / 1 2 



*\'<ci i/r+ / \z \'). (io) 

\JJQ JO 



Then, ([9]) is integrated over [0, T] and (|10|) is used to get 



> 2 <C[ II \f\ 2 + /Vol 2 ). Ill) 



Inequalities (|10[) and (|11|) finally imply that 



Q Jo 



\4 2 y <C[ I/ \f\ 2 + I No! 2 ). (12) 



IQ JO 

Equation Q is multiplied by (az xx ) xx and integrated over (0, 1). Some integrations by 
parts give also 



ld_ 

2 dt 



cr\z X x\ 2 )+J \{oz xx ) xx \ 2 



r 

I f(o~z xx *) xx 
Jo 
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We can then write 



Thus, we obtain 



d 



ilt y i cr\z xx \ 2 ^j + J <t\(z xx )xx\ 2 < J \f\ 2 - 



\\4 2 y 2 <c([[ i/r+ f L \ z ' ( ;\ 2 ). 

J Jo Jo 



(13) 



On the other hand, Equation (JS| is derived with respect to time. Thus q := z t satisfies 
qt + (<x {x)q xx )xx = ft, V{t,x)€Q, 



(14) 



«(t,0) = 0, 9(t,l)=0, Vte(0,T), 
qx(t,0) = 0, gb(t,l) = 0, VtG(0,T), 
_ g(0,x) = /(0,x)-(o-^'(x))'', Vie (0,1). 

Using estimate H3}, we obtain g G y 2 if (/(0, x) - (oz'o(x))") G # 2 (0, 1) and / t G 
L 2 (0,T;£ 2 (0, 1)). These hypotheses are fulfilled if z € H e f] H§(0,1) and / G T. Note 
that T C C([0, T]; H 2 (0, 1)). From the equation satisfied by z and the fact that / G T and 
zt G 3^2 , we determine that z G 3^6 , which concludes the proof of Proposition 12.11 □ 



Then, we focus on the linear problem with non-homogenous boundary conditions and 
low-order coefficients that depend on time. 

Proposition 2.2 Let z G H 6 (0, 1), f € T, G ly G 2 G // 1 (0, T; -ff 4 (0, 1)) and G # 2 (0,T) 
/or j = 1, ... ,4 satisfying the compatibility conditions with zo- Then, the equation 



(15) 



z t + (a{x)zxx)xx + ~/(x)zxx + G\z x + Giz = f, V(t, x) G Q, 

z(t,0) =hi(t), z(t,l) =h 2 (t), Vte(0,T), 

z*(t,0) = /i 3 (t), z*(t,l) =ft 4 (*)> Vte(0,T), 

z(0,a;) = zo(x), Vx e (0, 1), 

has a unique solution z G Z and there exists C > such that 



Proof. We first prove this result for null boundary data (i.e. for hj = for j = 1, . . . , 4 and 
therefore z e H 6 D H§(0, 1)). 

For any w £ Z, FL(w) is defined as the solution of (8} with / = (/ — 'y(x)w xx — GiW x — G2W). 
Note that / G J- and therefore n(tti) £ Z is well defined. 

If T is small enough, then II is a contraction. Indeed, for any w,w G 2, we have (the 
space L m (0,T; H n (0, 1)) is denoted as L m (H n )) 

\\H(w) - H(w)\\z < C\\y(x)(w X x - u>xx) + Gi(w m - w x ) + Gi(w - w)\\r 

< C\\W - 'W|| £ 2 (ir 6 ) + C\\Wt - Wt\\ L 2(H2) (16) 

< CT 3 \\W — »||i4( H 6) + CT* ||t« t — ™t||L 4 (H 2 ) 

< CT^IIlO-Mjllye+CrilK-tBtHys 

< CT*\\w-v)\\z. (17) 
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Hence, the operator II has a unique fixed point in Z, which is the solution of (|15[) with 
hj = for j — 1, ... ,4. Using standard arguments and the linearity of this equation, the 
solution can be extended to a larger time interval. 

In order to prove the general case, take hj £ H 2 (0,T), j — 1,...,4 compatible with 
2o. It is not difficult to find a function ip £ H 2 (0, T; C°°([0, 1])) satisfying the boundary 
conditions of (|15[1 . For instance take ip{x,i) = Y^j=i Pi ( x )hj(~t) where pi(x) = 2x 3 — 3a; 2 + 1, 
P2{x) = —2a; 3 + 3a; 2 , P3(x) = x 3 — 2x 2 + x and P4{x) = x s — x 2 . In particular we have 
Lip := ij)t + {o~(x)ip xx ) xx + 'y(x)?p xx + Gitp x + G21P £ T ' . Then, if w is the solution of equation 
(I15p with null boundary data, initial condition too — ?/>(•, 0), and right-hand side equal to 
/ — Lip, let us define z = w +ip. It is not difficult to see that z is the required solution. □ 

Remark 2.3 The third-order term z xxx can be added to Equation (|15f) . Indeed, in that case 
(|16p becomes C\\w — w \\l 2 (h 7 ) + C||wt — u)t\\ l 2 (h 3 )> which is bounded by 

CT* \\w - H|^i (H6) ||TO - ™lli2 2 (H 8) + CT z \\w t - w t ||^ (H2) ||w t - wt\\% 2 (Hiy 

This last expression is bounded by (|17p . The remainder of the proof is the same. 

Again, by using a fixed point theorem, we can prove Theorem 1 1.1 1 for equation ([7]). 
Let j/o £ H G (0, 1), hj £ H 2 (0, 1) compatible with yo, an d g £ J-. For any v £ Z, we define 
A(v) as the solution of ()15[) with / = (g — vv x ) and zq — yo- Note that / £ T and therefore 
A(t>) £ Z is well defined. Indeed, if v £ ^3 and t> t £ X)> then we have 

(vv x ') xxxx = (lOi'xx^'xxx ~t~ §v x v xxxx + f fxxx^x) £ L (0, Z/ (0, 1)) 

and 

{vv x ) t = v t v x + vv xt £ L 2 (0, T; L 2 (0, 1)). 
Furthermore, we can prove 

||A(«)|U < CI \\g\\ r + \\w w \\t + \\y \\ H e + ]T IIMh» 
{ i=i 

< CUgWjr + \\v\\% + \\y \\ He + £ \\hj\\ H 2 
{ 3=1 

Let e > to be chosen later and suppose that yo, hj and g satisfy @. Consider v such 

that \\v\\z < r with r > satisfying C(6e + r 2 ) < r. From (JTHJ) , we obtain ||A(ii)||^ < r. 

Thus, the application A maps the ball B r := {v £ Z/\\v\\z < r} into itself. 

We will now prove that A : B r — > B r is a contraction. For any z,v £ B r , A(z) — A(v) is 
the solution of (|15[) with 20 = 0, hj = for j = 1, . . . , 4 and / = uu^ — ^z^. We obtain the 
estimate 

\\A(z) - A(u)||.z < C\\vv x - zz x \\t < C\\(v - z)v x \\jr + \\z(v x - z x )\\t- 

Using the definition of the space F, v,z £ C([0, 1]; H 6 (0, 1)) ^ L°°(0, T; W 5 '°°(0, 1)) 
and v t , z t £ C([0, l];iT 6 (0, 1)) ^ L°°(0, T; W 1 ' 00 ^, 1)), we obtain 

\\A(z) - A(v)\\ z < C{\\v\\ z + \\z\\ z )\\v -z\\ z < 2Cr\\v - z\\ z , 
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which implies that A is a contraction if r is chosen small enough. More precisely, we can 
choose r, e such that 2CY < 1 and C(2e + r 2 ) < r. Hence, the map A has a unique fixed 
point y £ Z, which is the unique solution of (J7J). 
Thus, we have proven Theorem lf.il 



3 Global Carleman inequality 

In this section, a global Carleman inequality will be proven for the linearized K-S equation. 
We define the space 

V = {veL 2 {0,T;H 4 nH 2 (0, 1)) / Lv e L 2 ((0, T) x (0,1))} (19) 

where 

Lv — v t + {<TVxx)xx + qiVxx + qiVx + qov 
with qj £ L°°{Q) for j = 0, 1, 2. 

Consider /3 £ C 4 ([0, 1]) such that for some r > we have, for all x £ (0, 1): 

0<r<^(x), far* = 0,1, (20) 

§W<-<». (21) 

\<Tx{x)Px{x)\ < j mm .Mz)}. (22) 
4 ze[o,i] 



For instance, if a is constant, we can consider /3(x) = y/T+~x. 

On the other hand, given To G (0, T) we can choose <^o £ C 1 ([0, T]) such that 

9 i o (0) = ^o(T) = I and (23) 

< <t>a(t) < 0o (To) for each t G (0, T). (24) 

For example, if To = T/2, we can use 4>o(t) = t(T — i). 

We finally define the function 

4,(t,x) = £$r, (25) 
for (t, x) £ (0, T) x [0, 1], which is the weight function of the Carleman estimate and has a 
crutial role in the following result. 

Theorem 3.1 Let <f> be a function as \25\ and m > 0. Then there exists Ao > and a 
constant C = C(T, Ao, r, m) > such that if \\qi ||l°°((o,t)x(o,i)) ^ m f or * = 0, 1, 2 t/ien we 

-2A</> / u * + (CWicxjiCX . \7.7i |2 i \5;5i |2 . \3j.3i |2 . « ii i2\ j 

e v l — hA <p \v\ +\<p \Vx\ +A <p \v X x\ + \<p\Vxxx\ I dxdt 



JO 



T 



< C [ [ e-'^lLvf dxdt 
Jo Jo 

+ C j e- 2A0(t ' o) ^A 3 ^(f,O) ( 7(O) 2 |« a; 4f,O)| 2 +A^(t,O)a 2 (O)|^ M (t,O)| 2 ) dt (26) 
for all v £ V, /or aZi A > Ao . 
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Proof. Consider the following operator P defined in Wa := {e _Al *« : v G V} by 

Pw = e~ A ^L(e A, *w). 

We then obtain the decomposition Pw = P\w + P2W + Rw, where 

Piw = QX 2 4> 2 ow xx + X 4 </> x crw + (aw xx ) xx + QX 2 (4> 2 x a) x w x (27) 
P2W = w t + 4\ 3 4> x aw x + 4X(p x aw xxx + 4A 3 ' 4> x (4> x o) x w (28) 
Rw = \(f)tw + 2X(p x a xx w x + \ 2 (p x a xx w + X(f> xx a xx w 

+ 6\<j> x a x w xx + 6\ 2 (/> x (p xx cr x w + 6X(p xx a x w + 2X<p xxx a x w 

+ AX 2 (j} x (f) xxx cTw + 6X(/) xx aw xx + 3X 2 (/>l x aw + 4X(p xxx aw x 

+ X<f> xxxx ow + qow + qiw x + qiX<j> x w 

+ q-2W xx + 2Xq2(j) x w x + X 2 q 2 <t> 2 w + X</> xx q2W 

— 2A 3 4> 2 <j> xx aw — 2X i (f) x a x w. (29) 

Thus, 

\\Pw - Rw\\ 2 l2(q) = \\Piw\\ 2 L 2 {Q) +2{P 1 w,P 2 w) + \\P2w\\ 2 L 2 {Q) 

where (■, •) is the L 2 (Q) scalar product. 

For any v £ V we obtain v t £ L 2 (0,T;L 2 (0, 1)) and then v G C([0, T]; L 2 (0, 1)). From 
the construction of <j> (see (I23|> ). we obtain w G C([0, T]; L 2 (0, 1)) and w(x,0) = w(x,T) = 
for any w G W\. 

Let us define the notations 

I(w) = -6A 7 / / 4> x (j> xx a 2 \w\ 2 , 



I{w x ) = —X° J J <j) x tj(30<j) xx a + 12(j) x a x )\w x \ 
I(w2x) = —A 3 j (J (j}la(58(j> xx a + 4:0(j) x a x )\w xx 
I(w3x) = — A // o{2(j) xx a — 4:<f> x a x )\w xxx \ 2 . 



and 



f T 2 2 222 

h — (10A 4> x cr \w xx \ + 2X4> x aa xx \w xx \ +2X(f) x a \w xxx \ ) 

Jo 



The following weighted norm is defined, for any w G Wa, as 

\\w\\ 2 — \ \ (X 7 </)'\w\ 2 + X 5 <j) s \w x \ 2 + X 3 (j> 3 \w xx \ 2 + X(j>\w xxx \ 2 ) dxdt. 
Jo Jo 

We first require the following 
Lemma 3.2 Under the hypothesis of Theorem \3. 1\ there exists 8 > such that 

{P 1 w,P2w} L2(Q) >8\\w\\l^+I x (30) 

for X large enough and for all w G Wa ■ 
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Proof. It is sufficient to prove that 

3 

(Piw,P 2 w) L2 = ^2l(w kx ) +Ro(w) +I X 

for a large enough A, for all w £ W\, where |7io(w)| < A -1 ||w|| 2 . 

Let us first assume that we have pip. From hypotheses (|20[) to 
exists e > such that <f) satisfies for all x G (0, 1), 



(31) 



we know that there 



(32) 



4>xx(x) < —E(f) < 0, 

30<f> xx (x)a(x) + 12(/> x (x)a x (x) < -ecp < 0, 
58(p xx (x)a(x) + 40<f> x (x)a x (x) < —e<j) < 0, and 
2(f> xx (x)a(x) - 4</> x (x)cr x (x) < -e<f> < 0. 

Furthermore, from (|20p we can prove that <j> < C<j> x . Then from (|31[) we obtain, for a large 
enough A, 

3 

(P lW ,P 2 w) L2 =^2l(w kx ) + R (w) +I X 

> 28\\w\\ a x - \Ro(w)\ +I X (33) 

>*IMI» t+i*- 



Let us now prove (|31[) : we write {Piw,p2w) 



L2 (Q) = wnere denotes the 

L 2 -product between the i-th term of Pxw in (|27)l and the j-th term of P2W in (|28|l . 

Integrations by parts in time or space are performed on each expression Iij. Each 
resulting expression will be included in one of the terms on the right-hand side of ()31[) . The 
results are listed below, and we indicate for each term where it will be included. 



h. 



-h,i + 3A 2 // (4> x o) t \w x \ 2 



• 7i 2 = -12A J 



RoM 
i> x o 2 ) x \w x \ 2 . 



• h,3 



-12A J 



„cr ) x \w xx \ + 12A" 



j 3 2 1 1 1 



7i 4 = 12A J 



[<j> 3 x a((j)la) x ] xx \w\ 2 - 24A 5 // ^^(^cr)^!^^ 
Q J J Q 



• ^2,1 = 



• 7 2 ,2 



' 2 



-2A 7 [[ {<t>W) x \w\ 2 . 



I(w) 
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h.3 = -2A 



• 12,3 



• h,2 



• 7, 



5 II ^ cr2 ^ xxx ^ 2 + 6a5 jj (file 2 )* 



• 7 2 ,4 = 4A 7 



t> x a{cp x cr) x \w\ 



I(w) 



= 0. 



2A // [(^^cr)^^]^!™^! 2 +4A 3 



c a) x a\w x 



-f(™2x) 



+ 2A- 



c)iC0" |«k:a:| — 2A 



f 


J.3 2i i2 


/o 





2A 



f 


<t> X G<Tx X \w xx f' 




I 


/o 









^a xx ) x \w xx \ 2 + 8A ^o-crxliysxl^ 



+ 2A 



T 



2i i2 
C C7 U^mi 



^(™XXx) 

- 2A / / (ifcc^lwwl 2 - 
o •/ •'Q 



r{-w XX x) 



• 7 3 ,4 = 4A' 



jj (<t>x(4>l 



o) x ) xx aww xx — 4A 3 // (0 a; (<?f> 2 cr) I ) :I; CT|TO I | z 
Q 



Ko(«) 



H»(") 



+ 4A 



' jj 4>A<t> 2 x ^)x(y\w2 X \ z ■ 



• Ia,i = 6 A' 



ui t , which is canceled when adding with 7i, 



h,2 = 24A J / / (</> 2 cr) :z ^ 3 <7|i(j x | 2 . 



• 7 4 ,3 = 12A 3 



[(0x^)^0^^]x!c|^| 2 — 24A 3 jj (4>la) x 4> x a[ 



• 7 4 4 = -12A J 



/(™2x) 



Q 



R (w) 

Adding all the terms, we obtain ()31|) . 



□ 



Then, we will prove a Carleman inequality for the conjugated operator P. 
Lemma 3.3 There exists Xq > such that for all A > Ao we have, for all w £ Wa, 



10 JO 



(A 7 (/> 7 |u;j 2 + \°<j) J \w x \ 2 + A^^w^l 2 + A^lwxa^l 2 ) dxdt 

+ ll P lHli 2 (Q) + IIT^HIi^Q) < c\\Pw\\ 2 L 2 {Q) - I x 
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Proof. 

From hypothesis (|20[) and inequalities listed in (|32p . we know that there exists S > 
such that 

X)l(wto)>*||i<, (34) 

fc=0 

for a large enough parameter A. 

Besides, due to the definition of Rw in (|29|) and ||gi||L°o((o,T)x(o,i)) — m f° r £ = 0, 1,2 it 
is trivial to check that 



(35) 



II-RHIl2((o,t)x(o,u) < c ( a6 <^ 6 M 2 + ^ 2 ^ ^ 2 |tox| 2 + A 2 JJ^ <t?\w xx \ 2 

< CA^ 1 ll^ll 2 . 
Thus, for a large enough A, we have 

\\Piw\\ 2 l2 +2{P lW ,P 2 w) + \\P2w\\ 2 L 2 = \\Pw-Rwf L a 

< 2\\Pw\\ 2 L2 +2\\Rw\\ 2 L2 (36) 

< 2||P TO || 2 . 2 +CA- 1 |hll^. 

From Lemma [3. 2 1 and estimates (|36[) and (|34f) . we conclude the proof of Lemma 13.31 □ 
To complete the proof of Theorem l3.ll we have to deal with the norms fo Piw and P2W 
appearing in Lemma 13.31 From the definition of P2W, and because (|20[) holds, we have 

T-rlwil 2 < -T7 \P2w\ 2 + C (\°<f> a \w\ 2 + \°4> J \w x \ 2 + X(p\w xxx \ 2 ) 



and 



Jj) W A 2 < C II \P2W\ 2 + C\\w\\l^ 



IQ *V J J Q 

for a large enough A. A similar result is proven for (aw xx ) xx and Piw, and we then have 

-]-(\w t \ 2 + \{aw xx ) xx \ 2 ) <C [f (\Piw\ 2 + \P 2 w\ 2 ) +C|M| 2 10 - (37) 
From (I37|) and Lemma 13.31 we obtain 

^—(\w t \ 2 + \{uw xx ) xx \ 2 ) + \ 7 (f> 7 \w\ 2 + A 5 ^ 5 !^! 2 + \ 3 (f> 3 \w xx \ 2 + \4>\w xxx \ 2 dxdt 

Q X 9 

2 



<C J J \Pw\ 2 dxdt - CI X . (38) 

To handle the terms in I x , we note that for any x £ (0, 1) and a large enough A, 

-CX / (j) x (x,t)a(x)a xx (x)\w xx (x,t)\ 2 dt < CX 3 / tj) x (x, t) 3 a(x) 2 \w xx {x, t)\ 2 dt. 
Jo Jo 

Then 

r-T fT 

-CI X <CX 3 cj> x (0,t) 3 a(0) 2 \w xx (0,t)\ 2 dt + CX (j> x (Q,t)a{Q) 2 \w xxx {Q,t)\ 2 dt (39) 
Jo Jo 

and from (138 1) and (1391) we obtain 



— (\w t \ +\{ow xx ) 

XX \) + \M\\,<e<C 

Q X <P JJQ 

fT fT 

+ CX 3 / (t> x (0,t) 3 a(p) 2 \w xm (0,t)\ 2 dt + CX / cj) x {^t)a{0) 2 \w xxx {^t)\ 2 dt. (40) 
Jo Jo 
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Computing the derivatives of e^w it is trivial to prove that 



k 2 



= |d£(e A0 ™)| < C*^|A fc " j ft - J 9> 



for each fc = 0, . . . , 3. Therefore 

£ jT 1 e - 2A * (A 7 7 |e A ^| 2 + A 5 / Ke^w)*! 2 + A 3 3 |(e A0 w)«| 2 + A0| (e A *to)«„,| a ) dzdi 



Considering finally that Pw = e ^Lu, we obtain Carleman estimate (|26|l . □ 

Remark 3.4 VFe considered the function /3 to be increasing. This allows the Carleman 
inequality to be obtained with boundary terms at x = 0. If a decreasing function [3 was used 
instead, then an inequality with boundary terms at x = 1 would have been obtained. As 
discussed in the following section, the boundary terms in the Carleman inequality are related 
to the location of the observations in the inverse problem. 

4 Inverse Problem 

In this section, the local stability of the nonlinear inverse problem stated in Theorem 11.21 
will be proved following the ideas of [5] and [T3]. The proof is splited in several steps. 

Step 1. Local study of the inverse problem 

Let 7, 7, y and y be defined as in Theorem 11.21 If we set u = y — y and / = 7 — 7, then u 
solves the following K-S equation: 

u t + {<r{x)uxx)xx +-yu xx + yu x + y x u + uu x = f(x)y xx (x,t), V(t,x) G Q, 

u(t,o) = ti(t,i) = o, vte(o,r), 

(41) 

u w (t,0) =u x (t,l) = 0, Vte(0,T), 
u(0,x) = 0, Vz 6 (0, 1). 

Then, in order to prove the stability of the inverse problem mentioned in the introduction, 
it is sufficient to obtain an estimate of / in terms of u xx (-, 0), u xxx (-, 0) and u(Tq, •), where 
7 and y are given, 7 G U and it is the solution of Equation (|4ip . 

We begin by deriving Equation (|41[) with respect to time. Thus, v = ut satisfies the 
following equation: 



vt + (crv X x)xx + jvx X + yvx + y x v = fy xx t — g, V(i, x) G Q, 

v(t, 0) = v(t, 1) = 0, Vte(0,T), 

w !C (t,0)=« a! (t ) l) = 0, VtG(0,T), 

v(0,a:) = fR{x,0), Vase (0,1), 

where g(a:,£) = u(x,t)y xt (x,t) + u x (x,t)y t (x,t). 



(42) 
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The proof relies on the use of the Carleman estimate given in Theorem 13.11 This result 
will be used twice: First, Equation (|42[) allows to estimate v in terms of /, y xx and g; Then, 
Equation (|41|l will be used to handle the terms u and u„ which appear in the expression of 
the source term g. The details are given below. 

Step 2. First use of the Carleman estimate 

Similarly to the proof of the Carleman estimate, we set w = e _A< *u. Then, we work on the 
term 



1 = 2 I I w(t,x)wt(t,x) dtdx. 
Jo Jo 



10 Jo 

On the one hand, we can calculate / and bound it from below. Indeed, using w(0, x) 
'^°' x ^v(p, x) = for all x G (0, 1) and Equation ((41}, we can easily obtain 



I = 



f \w (T ,:r)| 2 dx 
Jo 



i 

e -2\^{T ,x) \{fy xx _ {au xx ) xx - JU XX - yu x - y x u - uu x ) (T ,x)\ 2 dx 



> / r^'^l/Wr |v«, (T„, ;r)| J dx-C\\u (T )||^ 4(0il) - C \\u (T )||* 1(0il) 



where C depends on ||7|| L oo(o,i), ||jf(To)|| wi,»(o,i) and ll°"llw 2 .°°(o,i)- 

On the other hand, in order to estimate I from above we apply the Carleman estimate (|40p 
to Equation (|42p using qo = y x and q\ = y, which are uniformly bounded in L°°((0, T) x 
(0, 1)) by the hypothesis in Theorem 11.21 We obtain that 
"1 r T a 



r* No 

1 = 21 / w(t,x)w t (t,x) dtdx 
Jo Jo 

< ((("A«M) W M)f*) 3 (£f ^gf 

< CA" 3 / /" e" 2A '*|/(a;)y :l;;ct (a;,f)| 2 dj;dt + CA" 3 f [ e'^^x,^ 2 dxdt 

Jo Jo Jo Jo 

fT 

+C\- Z / e" 2 ^ -*) (A 3 3 (0, t)a 2 (0) \v xx (0, t)\ 2 + X(f> x (0, t)a 2 (0)\v xxx (0, t) | 2 ) dt. 
Jo 

Step 3. Second use of the Carleman estimate 

Considering that g = uy x t+u x yt, we will now apply the Carleman estimate to Equation (|41[) 
in order to manage the term in g of the previous inequality. The unknown trajectory y is 
nevertheless such that y xt and y t belong to L°° (0, T; L°°(0, 1)) since y G H 1 (0,T; H 4 (0, 1)). 
Thus, we have 

// e- 2 ^\g(x,t)\ 2 dxdt < 2 [[ e- 2X *\u\ 2 \y xt \ 2 dxdt + 2 ff e - 2 ^\u x \ 2 \y t \ 2 dxdt 

JJQ JJQ JJQ 

< C [[ e" 2A0 (\u\ 2 + \u x \ 2 ) dxdt. 



The Carleman estimate (I40|l is applied to equation (|41|l . using the identity y x u + uu x = 
uy x , and taking go = y x and qi = y, which are bounded in L°°((0,T) x (0, 1)). 
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Therefore, we can choose Ao as large as possible in Theorem 13. II we then obtain 
'\g(x,t)\ 2 < CX~ 5 JJ e~ 2X * (\ 7 \u\ 2 + \ 5 \u x \ 2 ) 



e 



Q J J Q 

°- 2X *\f(x)y xx (x,t)f 

Q 

+ C\- 5 f e - 2A0(o ' t) (A 3 ^(O,t)a 2 (O)|^(O,t)| 2 + A0 :c (O,t)a 2 (O)|u M;c (O,t)| 2 )dt. 
Jo 

Gathering all of the estimates of / and g that were obtained above, we have 

f e- 2 ^ T '^\f(x)\ 2 \y xx (T , x)\ 2 dx-C\\u (T )|| 2 h4(0i1) - C \\u (T )\\% Hoj) 
Jo 

< C\- s JJ e- 2X4, \f{x)y xxt {x,t)\ 2 dxdt + C\- 6 JJ e- aA *|/(ar)y M (af,t)| a daf(ft 

j-T 

+ C\~ a I e- 2A0(o ' t) (A 3 ^(O,t)<J 2 (O)ii la;;c (O,t)j 2 + A0 ;c (O,t) ( 7 2 (O)|^ a; (O,O| 2 )^ 
+ CA -3 



/ e- 2A0(o ' t) (A 3 ^ 3 (O,t)<7 2 (O)!« ;ca; (O,t)| 2 + \M0, t)a 2 (0)\v xxx (0, t)\ 2 ) dt. 
Jo 

From the hypothesis of the theorem, we have y xx G L°°(0, T; W /1 ' oo (0, 1)), \y xx (T Q ,-)\ > r > 
in (0, 1). Using also that the Carleman weight function satisfies e' 2 ^^ < e - 2A ^( T o^) m 
(0,T) x (0, 1), we obtain 



e 

o 



2X ^ T ^\f(x)\ 2 dx 

< C(\- 3 J" e- 2 ^ T ^\f(x)\ 2 dx+\\u (T )|| 2 H4(0il) + ||u (T„)|| 4 fll(0il) 

+A- 8 f e ~ 2X< t ,i -°' t ) (A 3 3 (0, t)a 2 (0)\u xx (0, t)\ 2 + \<f> x (0, t)a 2 (0)\u xxx (0, t)\ 2 ) dt 
Jo 

JJ e- 2X ^°- t \\ 3 ^l(Q,t)a 2 (Q)\v xx {Q,t)\ 2 + A^(0, i)a 2 (0)1^(0, t)\ 2 )dt). 

Therefore, by choosing A large enough, we prove that there exists a constant C that 
depends on r, K, T, Ao, m such that 

\\m\\h (0 ,i) <c(n 

+ \\u xx (t, 0)||hi (0]T) + \\u xxx (t,0)\\ 2 H i i0iT) j. 

This estimate leads to the stability of the initial inverse problem and we have obtained 
the first estimate in Theorem 11.21 



Step 4. Regularity estimate 

We will now show that the second estimate also holds. Indeed, using To £ (0, T) and using 
that H^O.Tii/^O, 1)) ^ C{[0,T};H p (0, 1)), we have both, for some constant C, 

\\y(To,-) -y( T o,-)\\ H *( ,i) - G W y ~ yWHHo.T-.H^oA))' 

\\y(To, ■) - y(T , -)IIh1( ,i) ^ C\\y - y\\ H i (0yT . H i (0yl)) . 
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Since {y-y) € L 2 (0, T; H 4 n Hq(0, 1)), we have (y xxx - y xxx ) G L 2 (0,T; H 1 (0, 1)) and as 
L^TjiT^O, 1)) ^ L 2 ((0,T);C([0, 1])), we get the estimates 

\\y xx (-,0) - y xx (-,0)\\ L 2 (0 ^ T) < C||y - y|| L 2 (0jT;H 4 (0il)) , 

\\yxx X (-, 0) — y xxx (-, 0)|| L 2( ,T) < C||2/ — V\\l2(P,T;H*(P,1))- 

Using the same, but starting with (yt — yt) £ £ 2 (0, T; H 4 I1Hq(0, 1)), we obtain the estimates 

\\ytxx(;0) - ytxx(; 0)|U2 (0jT) < C||y t - y t || L 2 (0 , T;£r 4 (0il)) , 
||2/tx^(-,0) — ytxxx{-, 0)||i,2 (0)T ) < C\\yt — yt \\l*(p,t;h*(o,i)) ■ 

These inequalities provide the second estimate in ((6J and conclude the proof of Theorem 

mi. 
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